Magnification-Temperature Correlation: the Dark Side of ISW Measurements 
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Integrated Sachs- Wolfe (ISW) measurements, which involve cross-correlating the microwave back- 
ground anisotropies with the foreground large-scale structure (e.g. traced by galaxies/quasars), have 
proven to be an interesting probe of dark energy. We show that magnification bias, which is the 
inevitable modulation of the foreground number counts by gravitational lensing, alters both the 
scale dependence and amplitude of the observed ISW signal. This is true especially at high redshifts 
because (1) the intrinsic galaxy-temperature signal diminishes greatly back in the matter domi- 
nated era, (2) the lensing efficiency increases with redshift and (3) the number count slope generally 
steepens with redshift in a magnitude limited sample. At z > 2, the magnification-temperature 
correlation dominates over the intrinsic galaxy-temperature correlation and causes the observed 
ISW signal to increase with redshift, despite dark energy subdominance - a result of the fact that 
magnification probes structures all the way from the observer to the sources. Ignoring magnifica- 
tion bias therefore can lead to (significantly) erroneous conclusions about dark energy. While the 
lensing modulation opens up an interesting high z window for ISW measurements, high redshift 
measurements are not expected to add much new information to low redshift ones if dark energy is 
indeed the cosmological constant. This is because lensing introduces significant covariance across 
redshifts. The most compelling reasons for pursuing high redshift ISW measurements are to look 
for potential surprises such as early dark energy domination or signatures of modified gravity. We 
conclude with a discussion of existing measurements, the highest redshift of which is at the margin 
of being sensitive to the magnification effect. We also develop a formalism which might be of more 
general interest: to predict biases in estimating parameters when certain physical effects are ignored 
in interpreting observations. 

PACS numbers: 98.80.-k, 98.80.Es, 98.70.Vc, 95.36.+X, 98.62.Sb 



I. INTRODUCTION 

In a universe with accelerating expansion, the gravi- 
tational potentials associated with large-scale structures 
decay. A photon traveling through a decaying potential 
will experience a net change in energy. This leads to 
a secondary anisotropy in the cosmic microwave back- 
ground (CMB) called the integrated Sachs-Wolfe (ISW) 
effect [l|, [|[ . The ISW effect has been proposed as a tool 
to examine dark energy [3[ . Recent studies have detected 
the ISW signature providing an additional confirmation 
of the presence of dark energy 0, 0, i, 0, B H 0, [H O • 

Since inhomogeneities in the matter distribution trace 
gravitational potentials, one expects a significant correla- 
tion between the CMB temperature anisotropies and the 
distribution of large-scale structure at redshifts where 
there is a non-negligible fraction of dark energy. A uni- 
verse with non-negligible curvature or contributions from 
anything other than pressureless matter can do the same. 
Bearing in mind the evidence for spatial flatness [HI, [l4| , 
we will only consider dark energy. 

The ISW anisotropy can be isolated from the primor- 
dial anisotropies in the CMB by cross-correlating the 
temperature and galaxy/quasar fluctuations. (Hereafter, 
galaxy and quasar can be considered roughly synony- 



mous: essentially all statements about one apply to the 
other.) Cross-correlating the temperature anisotropy in 
direction 0, 5 T {6) = 6T(G)/T, with the galaxy fluctua- 
tion at redshift z in direction 6 , 5 n (6 , z) — 8n(0 , z)/n, 
gives 

(6 T (6)S n (6,z))=w nT {e,z) (I) 

where cos 8 = 60. This scale and redshift dependent 
signal provides information about the growth rates of 
large-scale structure at the redshift of the galaxy sample. 

Gravitational lensing alters the observed galaxy fluc- 
tuation in two ways. First, lensing alters the area of 
the patch of sky being observed thus modifying the ap- 
parent number density. Second, lensing can focus light 
promoting intrinsically faint objects above the magnitude 
threshold thus increasing the measured number density 
[l5l [lr3 |. Including these effects changes the observed 
galaxy fluctuation to 

S n = S g + <5 M (2) 

where S g is the intrinsic galaxy fluctuation and 8^ is the 
magnification bias correction due to gravitational lens- 
ing. The effect of magnification bias on the galaxy- 
galaxy and galaxy-quasar correlation functions is well 
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FIG. 1: The derivative (d/dz)[D(z)(l + z)] (solid line) where 
D is the linear growth factor, a galaxy selection function (dot- 
ted line) and the corresponding lensing weight function di- 
vided by c/Ho (dashed line). The galaxy-temperature cor- 
relation is proportional to the integral of the product of the 
solid and dotted lines, while the magnification-temperature 
correlation is proportional to the integral of the product of 
the solid and dashed lines. 

studied [M 03, E 0J, HI Ej]. The most recent mea- 
surements of this effect are discussed in [22, [52, [24| . Dis- 
cussions of earlier measurements can be found in the ref- 
erences therein. 

With magnification bias the cross-correlation signal 
becomes 

w nT {0, z) = w gT (0, z) +w^ T (0, z). (3) 

The galaxy-temperature term, w g T(0,z), is significant 
only when dark energy is non-negligible (assuming a flat 
universe, as is throughout this paper). Thus for high 
redshift galaxy samples this term is very small. However 
the magnification-temperature term, w^t{6,z)^ depends 
on the lower - redshift distribution of lensing objects so at 
high redshifts it may dominate over the w g T(9,z) term. 
Figure Q] illustrates this. 

We examine the effect of the magnification bias term on 
measurements of the ISW effect from cross-correlation. 
The magnification-temperature correlation has a dif- 
ferent scale and redshift dependence from the galaxy- 
temperature correlation. We will demonstrate that the 
magnification term can be large though its magnitude de- 
pends on the population of galaxies under consideration. 

This paper is organized as follows. We present the ba- 
sic expressions that govern the anisotropies and correla- 
tion functions in SjTj]and HID The ISW cross-correlation 
signal, especially in the presence of magnification bias, is 
sensitive to the sample of objects one is cross-correlating 



with. We define two illustrative samples in SjFV] The 
cross-correlation signals for these two samples, with and 
without magnification, are discussed in jJV] In WI[ we 
investigate how erroneous conclusions regarding dark en- 
ergy would be reached if one ignores magnification bias 
when interpreting ISW measurements. In Will we study 
how projections for the signal-to-noise and parameter 
estimation might be altered by the presence of magnifi- 
cation bias. We conclude with a discussion of existing 
measurements in Willi Appendix [A] is devoted to a 
technical discussion of how to estimate error and bias. 
In particular, we develop a formalism for predicting the 
estimation-bias for parameters of interest when certain 
physical effects (such as magnification bias) are ignored 
in interpreting data. We keep the discussion there quite 
general with a view towards possible applications other 
than ISW measurements. 

Two comments are in order before we proceed. Some 
of the existing ISW measurements come from cross- 
correlating the microwave background temperature with 
other diffuse backgrounds such as the X-ray background 
e.g. 0. Magnification bias, an effect that affects number 
counts, does not affect these measurements. This is not 
to say that lensing has no effect on diffuse backgrounds, 
gravitational lensing does have an effect through stochas- 
tic deflections. In particular, one might wonder how the 
lensing of the microwave background itself might impact 
the observed galaxy-temperature cross-correlation. The 
effect appears to be small on the large scales where the 
ISW measurements are generally considered interesting 
(we consider £ = 2 — fOO). The reader is referred to [J] 
and references therein for more discussions. 



II. GALAXY, TEMPERATURE AND 
MAGNIFICATION BIAS ANISOTROPIES 

The temperature anisotropy due to the ISW effect is 
expressed as an integral over conformal time from to 
today 770 [H| 

SJP W (6) = 2 r°drie-^^ (4) 
Jo or] 

where r(ry) is the optical depth between rj and 77. For 
perturbations sufficiently within the horizon, the gravi- 
tational potential <p is related to the mass (or matter) 
fluctuation 5 = 5p/p in Fourier space by (25| 

^(k,,) = -^O m (l + ,)^ (5) 

where il m is the ratio of the matter density to the critical 
density today, Hq is the Hubble constant today, c is the 
speed of light, z is the redshift, and k is the comoving 
wavenumber. The mass fluctuation S on sufficiently large 
scales grows according to linear theory: 6 oc D{z), where 
D(z) is commonly referred to as the growth factor. 

We are interested in cross-correlating the temperature 
anisotropies, St, with the observed galaxy overdensity, 
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S n . With magnification, the measured galaxy fluctuation 
is a sum of two terms (eq. [5] ) • The first is the intrinsic 
galaxy fluctuation 



(0) 



[6 M ] t (k,Zi) = 3Q ro ^(2.5 Si - 1) 



(12) 



6 g (6,Zi)= / dzb{z)W(z,Zi)6( X (z)e,z), 



where b{z) is an assumed scale-independent bias factor 
relating the galaxy overdensity to the mass overdensity 
i.e. S g — bS, W{z,Zi) is a normalized selection function 
about some mean redshift Zi, and x( z ) is t ne comoving 
distance to redshift z. 

The magnification bias term is also expressed as an 
integral over redshift [l6| 

6^(6, Zi ) = 30 m ^(2.5s(^)-l) 
c z 



X J dz J^j9(z,Zt){l + z)S(x(z)9,z). (7) 

The lensing weight function, g(z, Zi), can be thought of as 
being proportional to the probability of sources around 
Zi to be lensed by intervening matter at z [26| . 

g(z, Zi ) = X (z) rdz' x{zl) ~^ z) W{z', Zl ) (8) 
J z X[z ) 

For sources at a distance x, the lensing weight function 
peaks at ~ x/2. The prefactor of the magnification term 
depends on the slope of the number count of the source 
galaxies. This is defined as 



rflog 10 N(< to) 
dm 



(9) 



where to is the limiting magnitude of one's sample, and 
N(< m) represents the count of objects brighter than 
to. Magnification bias vanishes in the case that s = 0.4. 
The values of the slope s, as well as the galaxy bias b 
in eq. [6], depend on the population of galaxies under 
consideration. 



III. CROSS-CORRELATION FUNCTIONS 

For the most part, we focus on the spherical harmonic 
transform of the various angular correlation functions 
w §t{Q)i "UVt($) and so on. They can all be neatly de- 
scribed by the following expression: 

Cf B (zi,zj) = ^ J k 2 dkP{k)5Ai,{k 1 z, i )5B l {k,z 1 ) (10) 

where P(k) is the matter power spectrum today as a 
function of the wavenumber k, and the functions 8 At 
and SBg in the integrand can be one of the following: 

[S g ], (fc, Zi) =bi [dz W{z, z z )D(z) Je (k X (z)) (11) 



X ' dz Jj(^9(z,z t ){l + z)D{z) ]e {kx(z)) 



(ft) 



dz— [D(z)(l 

dz 



(13) 



z)] 



iijkxjz)) 
k 2 



The symbol ji denotes the spherical Bessel function. 
Here we have elected to assume a constant galaxy-bias 
bi and slope s, for each redshift bin centered around Zj. 
In the temperature integral we have neglected the factor 
containing the optical depth. This leads to an error of at 
most about 3% at the highest redshift we consider. 

Notice that the temperature function [5T]i(k) is in- 
dependent of redshift, so the cross-correlation functions 
Cf and Cf are functions of only a single redshift Zi- 
This is in contrast to the other correlation functions 
which depend on the redshifts of the two samples being 
correlated, e.g. C 9 t & '(zi, z 3 ), C^(Zi,Zj), C^ 9 (zi,Zj) and 
Cf fJ, (z i ,Zj). Note also that eq. [TU] and eq. [TJ] are not 
suitable for computing Cj T since that would account for 
only the ISW contribution to Cj T . 

Examining the correlation functions, we can see 
that the relative magnitude of the intrinsic galaxy- 
temperature correlation, Cf T , and the magnification- 
temperature correlation, C^ T , is redshift and scale de- 
pendent. The lensing efficiency increases with the red- 
shift of the source galaxies causing the magnification bias 
effect to increase with redshift as well. The scale depen- 
dence arises because each term depends on the matter 
distribution at different distances: the galaxy term de- 
pends on the matter distribution at the source redshifts 
while the magnification term depends on the distribution 
at the lenses. 

Additionally, the coefficients of C a g T and C^ T depend 
on properties of the galaxies via b and s. The magnitude 
of the magnification term compared with the galaxy term 
depends on the ratio (2.5s — Thus the correction 

may be negative for (s < 0.4) and in any case increases 
in magnitude as \s — 0.4 1 . The bias and slope depend on 
the choice of galaxy sample and are redshift dependent 
themselves. 

The Limber approximation, which is quite accurate 
when I is not too small (£ > 10), can be obtained from eq. 
[ID] by setting P(k) = P{k = (£ + l/2)/ X (z)) and using 
the fact that (2/tt) / k 2 dkj e (k X )ji(k X ') = (l/x 2 )S( X ~ 
X 1 ). We find that the substitution k = (£ + l/2)/x(z) 
is a better approximation to the exact expressions than 
k = £/x( z ) (see also @). More explicitly, Cf T and C% T 
under the Limber approximation are given by: 



n gT , \ 3fi m l?o bi 
[Zi) - 



c 2 (e + i/2y 



sH(z) , , 
dzW(z,Zi)—^D(z) 



d 



x — [D(z)(l + z)]P(k ± = (£+l/2)/ x ) (14) 
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Sample-I 














m < 27 


z-bin 1 


z-bin 2 


z-bin 3 


z-bin 4 


z-bin 5 


z-bin 6 


Zi 


0.49 


1.14 


1.93 


2.74 


3.54 


4.35 


k 


1.08 


1.37 


2.02 


2.90 


3.89 


4.81 


Si 


0.15 


0.20 


0.31 


0.43 


0.54 


0.63 


Hi 


34.6 


29.0 


10.1 


3.89 


1.68 


0.81 


Sample-II 














m < 25 


z-bin 1 


z-bin 2 


z-bin 3 


z-bin 4 


z-bin 5 


- 


Zi 


0.48 


1.07 


1.85 


2.67 


3.46 






1.13 


1.51 


2.73 


4.57 


6.63 




Si 


0.19 


0.35 


0.86 


1.31 


1.75 




Tli 


16.1 


8.6 


0.87 


0.11 


0.02 





TABLE I: The mean redshift z it the galaxy-bias bi, the slope 
Si and the number of galaxies per square arcminute rii are 
given for samples I and II, and for each corresponding redshift 
bin. 



Cf{z l ) = 



dzg(z,Zi)(l + z) 



x D(z) — [D(z)(l + z)]P(k ± = (t+ l/2)/ X ) (15) 

When displaying the above quantities in figures, we follow 
the custom of multiplying them by 2.725 K. 

Unless otherwise stated, we adopt throughout this pa- 
per the following values for the cosmological parame- 
ters when making predictions for the various correla- 
tions of interest: we assume a flat universe with a mat- 
ter density of fi TO = 0.27, a cosmological constant of 
Qa = 0.73, a Hubble constant of h = 0.7, a baryon den- 
sity of Slf, = 0.046, a scalar spectral index of n s = 0.95 
and a fluctuation amplitude of <jg = 0.8. The matter 
power spectrum is computed using the transfer function 
of Eisenstein and Hu [231 . The microwave background 
temperature power spectrum Cj T is computed using the 
publicly available code C AMB [28| , with an optical depth 
of t = 0.09 and no tensor modes. Throughout this pa- 
per, we consider ISW measurements at t = 2 — 100. In- 
cluding higher I modes does not significantly change our 
conclusions. This is because the signal-to-noise of ISW 
measurements is dominated by I = 10 — 100 (see e.g. 
[3ll l32j |. but also [HI). We adopt the Limber approxi- 
mation for all discussions in flVIl and IVIIl which concern 
issues related to signal-to-noise, estimation-bias and pa- 
rameter forecast, since these issues are not significantly 
affected by the very low I modes. The discussions and 
figures in EfVl however, make use of the exact expressions 
for the correlation functions (eq. [TU] h 



IV. SURVEY SAMPLES 

The predictions for the galaxy-temperature and 
magnification-temperature correlation depend on three 
quantities that are sample dependent: the redshift dis- 
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FIG. 2: The number of galaxies per unit redshift over the 
whole sky as a function of redshift, for two illustrative galaxy 
samples I and II (see text). The actual dN/dz for a given 
survey equals the sky coverage / s k y times the above. 



tribution of galaxies, the galaxy-bias, and the slope of 
the number count (at the magnitude cut-off). For illus- 
tration, we define two semi-realistic samples. 

Sample I is defined by an observed I-band (centered 
around 7994A) magnitude cut of 27 which implies a red- 
shift distribution dN/dz shown as a solid line in Figure 
[21 adopting the observed redshift-dependent luminosity 
function given by [29l |. This gives a net galaxy angular 
number density of 80 per square arcminute. We assume 
a sky coverage of f s k y = 0.5. These survey specifications 
are similar to those of the Large Synoptic Survey Tele- 
scope [30j . We divide these galaxies into 6 redshift bins, 
following the procedure of [311 ]: 



W(z,z{) oc 



1 dN{z) 

2 dz 



erfc 



erfc 



(i-l)A-z 
a(z)V2 
iA-z 



(16) 



(j{z)y/2 



where A = 0.8, a(z) — 0.02(1 + z), and Z{ denotes the 
mean redshift of the bin i. The complementary error 
function is defined as erfc(a;) = (2/y / 7r) J°° exp(— t 2 )dt. 
The normalization of W is fixed by demanding that 
J dzW(z, Zi) = 1. We have checked that none of our con- 
clusions are significantly altered by increasing the num- 
ber of bins. 

Given the observed luminosity function from [29], the 
slope s (eq. [5]) at the specified apparent magnitude 
cut can be obtained in a straightforward manner. The 
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FIG. 3: The galaxy-temperature correlation (solid lines) com- 
pared with the net galaxy-temperature plus magnification- 
temperature correlation (dotted lines) for four different red- 
shift bins. This is for sample I. 



galaxy-bias 6, on the other hand, requires some theoret- 
ical input. Here, we adopt the method of Kravtsov et 
al. [3J|. It works basically by number matching: one 
matches the the number of galaxies brighter than some 
magnitude with the number of halos and subhalos above 
some mass threshold. Once this mass threshold is known, 
the galaxy-bias can be computed. More specifically, the 
number of (central plus satellite) galaxies within a halo 
of a given mass M is 



■ erfc 



lnAfn 



M 



O.Ia/2 



lnM\ , 

) x(1 + 2om; 



) 



(17) 



Here, the complementary error function gives essentially 
a step function with a modest spread (of 0.1), and M can 
be thought of as the mass threshold. The quantity 20Mo 
specifies the (parent) halo mass at which the expected 
number of satellites is unity. This factor is in princi- 
ple redshift dependent, but we ignore such complications 
here. We determine M by demanding that the integral 
of the halo occupation number over the halo mass func- 
tion equals the number of galaxies for the sample under 
consideration. The bias of this sample of galaxies can 
be obtained by integrating the product of the halo oc- 
cupation number and the halo-bias over the halo mass 
function, normalized by the total number of galaxies. We 
adopt the mass function and halo-bias of Sheth and Tor- 
men [361 ] . 

Sample II is defined in a similar way to sample I, except 
that the limiting I-magnitude is taken to be 25 instead of 



FIG. 4: The analog of Figure [3] for sample II. 



27. The corresponding dN/dz is shown as a dotted line 
in Figure [2] The net galaxy angular number density is 26 
per square arcminute. We assume the same sky coverage 
of fsky = 0.5. This is a shallower survey, and we divide 
the sample into only 5 redshift bins i.e. apply eq. |16| up 
to i = 5. 

A summary of the mean redshift, the galaxy-bias b, the 
slope s and the galaxy number density for each redshift 
bin and each sample is given in Table HI The tabulated 
dN/dz for each sample is available from (35j . 



V. THE CROSS-CORRELATION SIGNAL 

In Figures [3] and 2] we show the cross-correlations 
Cf T (zi) and C" T (zi) = C 9 ( T [zi) + C^ 7 \zi) at several 
redshift bins for both sample I and sample II, calculated 
using equations (fl"0)) - (fl"3j) . 

There are several features in Figures [3] and 21 that 
are worth commenting. As is expected, the galaxy- 
temperature correlation Cf T decreases rapidly with red- 
shift: at high z, dark energy becomes subdominant mak- 
ing the time derivative of the gravitational potential quite 
small (eq. [4]). Note that this decrease of C| T occurs de- 
spite the overall increase of the galaxy-bias with redshift. 



The magnification-temperature correlation Cg T , on 
the other hand, generally increases with redshift be- 
cause of the increase of the lensing efficiency. The 
magnification-temperature correlation C mT can be high 
even as C gT becomes small, because the rather broad 
lensing weight function makes C^ T sensitive to structures 
at low redshifts where the gravitational potential has a 
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FIG. 5: The 2-point galaxy-temperature correlation 
w g T(0,Zi) (solid lines) compared with the net galaxy- 
temperature plus magnification-temperature correlation 
w n T(Q,Zi) (dotted lines) in several of our redshift bins. This 
is for sample I. 



significant time derivative. The general increase of Cf 
with redshift is aided also by the increase of the number 
count slope s — at a higher z one is generally looking at 
intrinsically brighter galaxies which reside in the steep 
part of the luminosity function. The net result is that 
the total cross-correlation C™ T starts to climb with red- 
shift, for z > 2. In fact at a sufficiently high z, C™ T 

even becomes comparable with C™ T or Cf T at the low- 
est redshift bin. Comparing sample I and sample II, one 
can see that a brighter magnitude cut leads to a more 
pronounced magnification bias effect. This is due to the 
steeper number count slope. 

Another feature visible in Figures [3] and Q] is that mag- 
nification changes the shape of C" T with I, The magni- 
fication term, being sensitive to structures at lower red- 
shift, peaks at a lower I than the galaxy term. For posi- 
tive magnification contribution, the peak of C™ T is at a 
lower I than the peak of Cf T . 

In Figures [5] and [6] we plot the real space correlation 
functions w 9 t(0, Zi) and w n T(0, Zi) without the monopole 
and dipole contributions. These are calculated from 



w x (0,Zi) 



500 

E 

1=2 



2£ + l 

47T 



Cf(z t )P £ (cos6). (18) 



where X symbolizes gT or nT and Pi(cos9) are the Leg- 
endre polynomials. We have extended the range of the 
sum to £ = 500 to insure the convergence of the sum for 
small values of 6. The general discussion about C^ T ap- 
plies to w n T as well; at z > 2 magnification bias grows to 
become an important component of the cross-correlation 
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FIG. 6: The analog of Figure [5] for sample II. 



signal, and the amplitude is preferentially increased at 
large angular scales. 

Magnification bias can of course have either sign. At 
low redshifts, one is generally looking at intrinsically faint 
galaxies which leads to a small s that can be less than 0.4. 
This is why Cf T < Cf T at low redshifts for our samples. 
The transition redshift from s < 0.4 to s > 0.4 is sam- 
ple dependent. For a sample with an apparent limiting 
magnitude cut, the brighter the cut, the lower the tran- 
sition redshift. It is conceivable that a sufficiently faint 
cut can be achieved such that s remains less than 0.4 
out to redshifts where \C^ T \ > C^ T . In such a case the 
net cross-correlation C™ T would become negative. This, 
however, requires an exceptionally deep survey: an I- 
band magnitude limit of 28 or higher, according to the 
observed luminosity function [29| . Conversely, if a suffi- 
ciently bright galaxy cut were used at low redshift it is 
possible that magnification could become important at 
lower redshifts than evidenced by our samples. 

In the literature the cross-correlation signal is of- 
ten presented in the normalized form Cf T (zi) / y/Cf 9 (zi) 
or w q t(9, zi) I \Jw gg (0, zi) to remove the dependence on 
galaxy bias. However, when magnification bias is taken 
into account the auto-correlation is modified to be 

Cp(z i) z j )= (19) 
Cf 9 (zi 1 Zj) + (zi, Zj) + ' (z{, Zj) + (Jf 1 (zj . Zj). 

Note, then, that the galaxy-bias would be completely re- 
moved by the division only if magnification bias were ab- 
sent. Had we plotted the normalized versions of the cross- 
correlation signals, the fractional difference between 
Cf Is/Cf and C\ lT ' l^fC™ would remain about the 
same as that between the un-normalized versions. This 
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FIG. 7: The average inferred dark energy equation of state 
w from ISW measurements when (erroneously) ignoring mag- 
nification. The true w equals —I. The solid symbols show 
the inferred w for each redshift bin separately. The dotted 
symbols at z = show the inferred w combining all redshift 
bins. The error-bars are la — note that they are dependent 
upon the inferred model (see text). The upper panel is for 
sample I, and the lower panel is for sample II. 



is because for the cases considered here C™{zi 
Cjl 9 (zi,Zi) to within about 20%. 



VI. ESTIMATION-BIAS DUE TO IGNORING 
MAGNIFICATION — - A THOUGHT 
EXPERIMENT 

Clearly, magnification bias should be taken into ac- 
count when interpreting ISW measurements, especially 
at high redshifts where the magnification-temperature 
correlation actually dominates over the usual galaxy- 
temperature signal. Ignoring it would lead to erroneous 
conclusions about dark energy. We quantify this by com- 
puting the estimation-bias in the dark energy equation of 
state w if magnification were ignored when interpreting 
measurements from samples I and II. 

A thought experiment is set up as follows. Suppose 
the true dark energy equation of state is w = — 1 i.e. a 
cosmological constant (and other cosmological parame- 
ters take the values stated at the end of ^IIII) . Suppose 
further one were to infer w (and other parameters) from 
measurements of the total cross-correlation C™ T by fit- 
ting them with a model that ignores magnification i.e. 
Cf T , and one were to assign errorbars based on the (er- 
roneous) model with no magnification. The fitted pa- 
rameters include: the dark energy equation of state w 
(i.e. w = P j p where P and p are the dark energy pres- 



sure and energy density respectively) , the matter density 
Q m (we assume a flat universe so that the dark energy 
density is 1 — Cl m ), the amplitude of fluctuations as, the 
spectral index n s , the Hubble constant h, the baryon 
density fif, and the galaxy-bias for each redshift bin bi. 
Constraints from other data (such as from the temper- 
ature anisotropy itself Cj T and from the galaxy power 
spectrum) are simulated by the inclusion of priors: the 
fractional (1 a) errors on f2 m /i 2 , flbh 2 , h and as are as- 
sumed to be 5%, the fractional error on bi is 10%, and 
n s has an (absolute) error of 0.02. These are similar to 
current constraints, depending somewhat on assumptions 
[l3T ] . A simple analytic expression can be derived for the 
parameter estimation-bias, which is discussed in detail 
in Appendix A (eq. }A14j ). 

The inferred w, together with its marginalized error, 
for samples I and II are shown in Figure [7] The modes 
used range from I = 2 to I = 100. The solid symbols 
show the inferred w from each redshift bin separately. 
The sign of the resulting estimation-bias is determined 
by the slope s. If s < 0.4, as is the case at low redshifts, 
the inferred w tends to be lower than the true value. 
This is understandable because the observed C nT is sup- 
pressed by magnification in that case, and the only way 
to stick with an erroneous model without magnification is 
to resort to a lower w to delay dark energy domination. 
The reverse happens at high redshifts, where s > 0.4. 
In fact, at the highest redshift bin, the bias can become 
quite large. Ignoring magnification bias is simply unac- 
ceptable. 

It is worth emphasizing that the errorbars are model 
dependent i.e. they are computed using the inferred 
model, which is different for each redshift bin. The gen- 
eral trend is for the errorbar to shrink as one goes to 
higher redshifts, despite the drop in signal-to-noise (see 
flVH)) . This is due to the fact that dhiC 9T /dw rises with 
both z and w. The generally positive bias in w at high 
redshifts helps to diminish the corresponding errorbars. 

The dotted symbols at z = show the inferred w and 
its associated errorbar when measurements from all the 
redshift bins are combined. The previous remarks on the 
model dependent nature of the errorbar applies here as 
well. In the case of sample II, one finds that the overall 
inferred w is biased high by almost 3 a. 



VII. SIGNAL-TO-NOISE AND PARAMETER 
FORECAST 



The above considerations suggest that magnification 
could cause the signal-to-noise of ISW measurements to 
remain favorable even at relatively high redshifts, deep 
into the matter dominated regime. This is borne out by 
Figure [5] [13] , which shows (with / s k y divided out) 
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FIG. 8: The signal-to-noise (squared) of ISW measurements 
from each redshift bin for samples I and II. The solid line 
ignores magnification (equivalent to setting s = 0.4) while 
the dotted line includes magnification. Note (S/N) 2 oc / sky 
i.e. the division by / s k y takes out the dependence on sky 
coverage. 



FIG. 9: The net cumulative signal-to-noise (squared) of ISW 
measurements from all redshift bins up to z max , for samples 
I and II. The solid line ignores magnification (equivalent to 
setting s = 0.4) while the dotted line includes magnification 
(the two are very close to each other). The division by / s k y 
removes the dependence on sky coverage. 



where rn is the mean surface density of galaxies, and we 
include modes from I — 2 to £ — 100 (as in the rest of the 
paper). The result for sample II is particularly striking: 
the signal-to-noise remains more or less flat out to high 
redshifts (dotted line), in sharp contrast with the signal- 
to-noise if magnification is ignored (solid line; equivalent 
to putting s = 0.4, or replacing the superscript n by g in 
cq. [20]). 

Magnification bias, therefore, opens up a high redshift 
window for ISW measurements which otherwise would 
not exist. That is the good news. The bad news, how- 
ever, is that since the magnification of redshift Z{ galaxies 
probes structures at redshifts z < Zi, the C™ T measure- 
ments at high redshifts are in fact quite correlated with 
those at low redshifts. Taking into account such correla- 
tions, we show in Figurc[9]thc net accumulated signal-to- 
noise (squared) for measurements from all redshift bins 
up to a given z max , defined as follows: 



N 
where 



(^max) 



2,<2 m „ £ 

(21) 



[Cavi]^ = 

(Cpjzj, Zj ) + Sq/njef T + Gf^Cf^) 

/sky(2* + l) 



(22) 



It appears the cumulative signal-to-noise reaches a 
plateau by z ~ 2 whether or not magnification bias is 
included. As is expected, magnification bias makes much 
less of a difference in this case: the solid and dotted lines 
are very close together. This happens despite the signif- 
icant difference in signal-to-noise on a redshift bin by 
redshift bin basis, shown in Figure [8] The culprit is the 
significant correlation between different redshift bins in- 
troduced by lensing — the high redshift measurements 
are in fact only probing large-scale structure at low red- 
shifts. 

The reader might also wonder why the two samples 
have such similar cumulative signal-to-noise. This is 
because at low redshifts (where most of the cumulative 
signal-to-noise comes from), neither sample is shot- noise 
dominated i.e. for the most part, the noise is dominated 
by the term associated with C™ n Cj T (eq. [2"2"]). 

Ultimately, we are interested in cosmological con- 
straints from ISW measurements on, for instance, the 
dark energy equation of state w. It is therefore useful to 
show the expected errors for w, both on a redshift bin by 
redshift bin basis (as in FigureJS) and on a cumulative 
basis (as in Figure (see also |33|). This is displayed in 
Figures [TO] and EJ 

The fiducial cosmological model here, as is in the case 
of the computation of S/N, is that described at the end of 
mill The cosmological parameters that are marginalized 
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FIG. 10: The marginalized error on w, ^(5w 2 ), for ISW 
measurements from each redshift bin. The solid line ignores 
magnification (equivalent to setting s = 0.4) while the dotted 
line includes magnification. 



FIG. 11: The marginalized error on w, ^{Sw 2 ), for ISW 
measurements from all redshift bins combined up to « mw . 
The solid line ignores magnification (equivalent to setting s = 
0.4) while the dotted line includes magnification (the two are 
very close together). 



VIII. DISCUSSION 



over and the set of priors assumed are identical to those 
described in §VII (with the addition of Si, the number 
count slope for each redshift bin, as a new parameter, for 
which we impose a prior of Asi — 0.02), and the Fisher 
matrix formalism for doing so is explained in Appendix 
A, in particular eq. |A17j . Once again, we see that while 
the presence of magnification bias improves the errors at 
high redshifts on a redshift bin by redshift bin basis, it 
does not in fact add much cumulative information if one 
considers the errors from combining redshift bins. We 
have experimented with changing the priors and varying 
sample definitions, and it appears this conclusion is quite 
robust. 

The most compelling reason for doing high redshift 
(z > 2) ISW measurements is therefore not that the 
expected cumulative error on w would continue to im- 
prove, but that there might be surprises e.g. dark energy 
might be quite different from the cosmological constant 
and actually remain a significant component of the uni- 
verse even at high redshifts, or gravity ma y be modified in 
non-trivial ways on large scales [38|, [39|, 14CJ, |4l|, |44| . At the 
very least, high redshift ISW measurements constitute a 
consistency check that one should make. As explained 
in £|VI1 such high redshift measurements should be inter- 
preted with care: magnification bias must be taken into 
account. 



A natural question: to what extent should we worry 
about magnification bias when interpreting existing ISW 
measurements? This is partially addressed in Figure [TSJ 
which shows the cross-correlation signal at 6° separation. 

The points with errorbars are measurements taken 
from a compilation in 42] with the addition of recent 
results from [ijj and [nf. A11 involve correlating the 
microwave background with galaxies, except the highest 
redshift one [llj which involves correlation with quasars. 
The median redshifts of the measurements are: z ~ 
0.1,0.15,0.3,0.5,0.9,1.5. The ISW signal shown is nor- 
malized by the galaxy/quasar bias, with the bias values 
taken from [13, EE SI - They are, in the order of increas- 
ing redshift: 1.1, 1.0, 1.2, 1.4, 1.7, 2.3. 

The solid lines show the theoretical prediction for 
Wgr/b (i.e. ignoring magnification) adopting the cosmo- 
logical parameters stated in mill Because the selection 
function is somewhat uncertain, we choose to illustrate 
the prediction with a 'broad' and a 'narrow' selection 
function, parametrized as follows: 



W{z) 



r 



q+l 



(23) 



where a = 2, /3 = 1.5 is broad, and a — 2 (1+3 zq), [3 = a 
is narrow. For each selection function zq is adjusted to 
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FIG. 12: The two panels show the 2-point correlation func- 
tion divided by the galaxy bias, w(0,z)/b(z), at 9 — 6° as a 
function of the median redshift of the sources, z. The points 
with errorbars are the existing measurements. The solid line 
shows the theoretical prediction ignoring magnification for a 
flat cosmological constant model (see £]III[) . The dotted lines 
include magnification for the cases s = 0.8 (upper dotted line) 
and s — 0.2 (lower dotted line). The upper panel uses a broad 
galaxy selection function while the lower panel uses a narrow 
selection function (see text). 



reproduce the median redshifts of the measurements. 

The dotted lines [43[ in Figure [T2l show the theoretical 
prediction for w n T/b (i.e. accounting for magnification) 
for s = 0.8 (upper) and s = 0.2 (lower). These values for 
the number count slope seem to span the range observed 
in quasar samples [23[ similar to the one used in the high- 
est redshift ISW measurement, where magnification bias 
is most relevant. 

Considering all angular scales and both selection func- 
tions, we find that for the four lowest redshift measure- 
ments \Wfj,T /w 9 t\ < 0.1, unless s > 1.4 or s < —0.6. For 
the NVSS sample at z ~ 0.9, we estimate the slope of 
the NVSS radio sample to be s = 0.320 ± 0.07 which 
results in Iiw^t/uVtI < 0.1. This is small compared 
with the measurement error. It appears that the quasar- 
temperature correlation at z ~ 1.5 is the only measure- 
ment for which magnification could be significant. To 
determine precisely to what extent this measurement is 
affected by magnification requires a more detailed analy- 
sis of the quasar sample. It is conceivable that a suitably 
chosen quasar subsample could exhibit a large magnifi- 
cation effect. 

Figure[l2]also shows that for z > 1 there is a significant 
dependence on the width of the selection function. This 
is because the mass-temperature correlation decreases 
rapidly with redshift. If a broad selection function is 



used, the galaxy-temperature correlation will decline less 
rapidly because low redshift sources will contribute. On 
the other hand, w^t depends on the lensing weight func- 
tion (eq. [5]) which is broadly distributed regardless of 
the width of the source distribution. In other words, 
changing the width of the selection function from broad 
to narrow causes the galaxy-temperature term to de- 
crease, while leaving the magnification-temperature term 
largely unchanged. Figure [12] illustrates this for a fixed 
9, however, the statement holds for other angles. Source 
distributions used in the existing ISW measurements are 
probably closer to the broad selection function. 

To conclude, we find that magnification bias alters the 
observed galaxy/quasar-temperature cross-correlation 
significantly at high redshifts (Figures [3] and [3}. The 
precise magnitude of the modification depends on the 
sample of objects under consideration. Three facts more 
or less guarantee magnification plays a non-negligible or 
even dominant role at z > 2: the generally steepening 
slope of the number count at high redshifts (because one 
is looking at intrinsically brighter objects), the decay of 
the intrinsic galaxy-temperature correlation C gT as one 
enters the matter-dominated era, and the increase of the 
lensing efficiency with redshift. Ignoring magnification 
bias when interpreting high redshift ISW measurements 
would lead to erroneous conclusions about the nature 
of dark energy. For instance, the estimated equation of 
state w can differ from the true one by more than 3 a in 
some cases (Figure [7]). 

The boosting of the ISW signal by magnification (as- 
suming s > 0.4 which is generally the case at high red- 
shifts when one is looking at intrinsically bright objects) 
implies that, despite naive expectations, ISW measure- 
ments remain viable even at z > 2. However, because 
of the correlated nature of the lensing signal (across dif- 
ferent redshifts), the cumulative information from low 
z's to high z's is not significantly enhanced by magni- 
fication. The most compelling scientific justification for 
pursuing high redshift ISW measurements is to look for 
surprises: for instance, dark energy might remain signif- 
icant out to high redshifts. At the very least, one would 
like to perform a consistency check of the very successful 
cosmological constant model of dark energy. The large- 
scale nature of the ISW signal also means it is a natural 
place to look for spatial fluctuations in dark energy [3l[ 
and possible signatures of modified gravity [44j . 

When interpreting future high redshift ISW measure- 
ments, accounting for the effect of magnification bias is 
a must. The net cross-correlation signal depends not 
only on the cosmological parameters (e.g. O ro , £Ia, h, 
n s , as) but also on the sample-dependent parameters b 
and s. The slope s can be estimated from the observed 
number counts. The bias b can be inferred by comparing 
the amplitude of the galaxy/quasar auto-correlation with 
the amplitude of the mass auto-correlation predicted by 
the microwave background, modulo the fact that the ob- 
served galaxy/quasar auto-correlation function is itself 
affected by magnification bias (though it appears to be at 
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most a 20% effect at the scales and redshifts we study). 
Alternatively, the bias can be estimated by examining 
the galaxy/quasar higher order correlations. The up- 
shot is that separating the galaxy-temperature and the 
magnification-temperature signals is in principle possi- 
ble, at least in a model dependent manner. Other possi- 
bilities are to exploit the sample-dependent nature of b 
and s to study the variation of the total ISW signal with 
different sample cuts, or to use the galaxy-galaxy correla- 
tion across different redshift bins to check for consistency 
with the ISW-magnification bias signal. 

An interesting question remains to be explored: given 
a survey with many different kinds of objects, such as 
galaxies/quasars of different luminosities and types, how 
should one go about choosing a weighting scheme that 
optimizes the total signal and/or its magnification com- 
ponent? The magnification component might be partic- 
ularly interesting if the galaxy/quasar bias turns out to 
be uncertain. We leave this for future investigation. 



Assuming Gaussian distributed data, the likelihood is 
proportional to exp [— x 2 /2], where 



X 



= Y J {di-d i )Cr j \d j -d j ) 



(Al) 



where dj = ((di — di)(dj — dj)), with di = (di). The vec- 
tor (component) di is in general some complicated func- 
tion of the parameters p a . The trick is to turn this into 
a linear estimation problem by supposing that p a is close 
to some fiducial value p a , and therefore: 



di = di + ^-5p 



Sp a =Pa~Pa 



(A2) 



where di is the expected data average if the parameter 
values were indeed p a , and the partial derivative is eval- 
uated at p a = p a - 

Substituting the linear model eq. |A2| into the expres- 
sion for x 2 > we obtain 
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APPENDIX A: ON ESTIMATION BIAS AND 
ERROR 

The technology for making error forecasts is well doc- 
umented in the literature [25j . What seems to be less 
often discussed is the issue of predicting the bias: what 
is the resulting bias in the estimation of a parameter if 
some particular physical effect is ignored or some ap- 
proximation is made in modeling the data? In our case, 
the parameter of interest is w and the physical effect 
is the magnification bias. We keep our discussion here 
fairly general since it might be of interest to researchers 
working on other problems. Much of our discussion is 
a straightforward generalization of the classic paper by 
Rybicki and Press [461 ]. 

In general, we have the problem of trying to obtain con- 
straints on m parameters, labeled p a with a = 1, 2, m, 
from an n-dimensional data vector, di with i = 1,2, n. 



X 2 = }_^(Spa - SpmFa P (6 P0 - (A3) 

a,/3 

up to an additive constant independent of Sp a - Here, 

SpT x = £ Fj^-^id, - dj) (A4) 

where the " on top of 5p™ ax is there to remind us that 
this quantity depends on the data dj (i.e. a stochastic 
quantity), and F a p is the Fisher matrix 



_ x ^ ddi t ddj 

1,3 



(A5) 



Therefore, the maximum likelihood estimate for p a is 

(A6) 



pT x =P a + 5Pa ! 



and the error covariance associated with such an estimate 
is given by the inverse of the Fisher matrix: 

(ffT) ~ (Pa (pT^) = F J (A7) 

Combining eq. |A4j and |A6| . we see that 

(Pa aX > = Pa + E W °i(.&) ~ ( A§ ) 
3 



04 



max\ ^.true 



Pa 



The estimation bias of interest is therefore (p; 
where Pa UC is the true value of the underlying parame- 
ter. Note that both the weights w a j and the vector dj 
are model dependent. For instance, one could choose to 
ignore certain physical effects, such as magnification bias, 
in constructing these quantities. The result would be a 
biased estimate of p a - 
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Let us suppose first that one is using the correct 
or exact model. Then, substituting (dj) = dj + 
[ddj/dp^ip*™ -p 7 ) into eq. [58], we have (pT*) = 
p^ rue , as expected. 

Suppose instead one is using an incorrect or approxi- 
mate model, where a certain physical effect (e.g. mag- 
nification bias) has been ignored. This means dj and 
the weights w a j are all computed by ignoring this ef- 
fect, whereas (dj), which is the expectation value of the 
observed data, of course has such a physical effect (and 
all other relevant physical effects) taken into account. 
The preceding argument (concerning the exact model) 
no longer works and one is left with an estimation bias 
i.e. (p™ ax ) 7^ Pa UC - Note that this bias depends on the 
assumed fiducial parameter value p a (as opposed to the 
case of unbiased estimation, where (p„ ax ) does not de- 
pend on p a at all). Generally, a modeler would choose 
a fiducial value that is close to the perceived (i.e. possi- 
bly incorrect, due to the use of an incorrect/approximate 
model) maximum likelihood value. We therefore obtain 
the estimation bias by iteration: starting with some fidu- 
cial parameter value, compute (p™ ax ), then adopt the 
result as a new fiducial value and iterate. In all cases 
investigated in this paper, a few iterations are sufficient 
to guarantee convergence. The modeler would assign a 
(theoretical) errorbar to his/her biased estimate based on 
eq. |A7j . where the Fisher matrix is, again, that of the 
approximate model. 

Lastly, let us address the question of priors. Priors can 
be thought of as just another kind of data, and they add 
to the x 2 m eq. |A1| a term of the form: 

E (ft* - Pa) M a0 (P0 ~ Pff ) 

a, 

where p a denotes the assumed prior value (e.g. fl m = 
0.27, etc) and M a p is the covariance associated with this 
set of priors. Completing square as before, one obtains, 
up to an irrelevant additive constant, a \ 2 slightly differ- 
ent from the one in eq. |A3j : 

X 2 = $>p„ - 6p™*)F^(6 Pf3 - (A9) 

a,0 

where F' a p is a generalized Fisher matrix, which is related 
to the original Fisher matrix by 

F^ = F a p + M-£ (A10) 

and the maximum likelihood estimate p™ ax — p a + <5p" lax 
is given by (superseding eq. | A4j ) : 

S0T - E F 'al [E - ^j) (All) 

+ E M ^7-P7) 

7 

Therefore, the expectation value of p™ ax is given by 

(c ax > = p a + E F '«l [E f^ 1 ^) - ^)( A12 ) 

id Pl3 



+E ik W((P7> -Pi) 

7 

A set of priors that are unbiased implies (p 7 ) = p^ rue . 
The estimation bias is given by (p'° ax ) — Pa 110 , an d the 
error covariance is given by the inverse of the generalized 
Fisher matrix: 

/"max ~max\ / -max\ / --max\ t^i~ 1 /aio\ 
(Pa P0 )-\Pa ) VP/3 ) = F a0 ( A13 ) 

Eq. |A12j and |A13| constitute the main results of this 
Appendix. Let us illustrate their use by applying them to 
a question of interest in this paper: what is the resulting 
bias on cosmological parameters if one ignores magnifi- 
cation bias when interpreting ISW measurements? The 
parameters p a stand for w, Q m and so on. The data 
(dj) are the C gT + C mT at various £'s and redshifts (col- 
lectively labeled by the index j here). The Cy matrix 
here is the covariance described in eq. [5^]. To figure 
out the estimation bias on parameter p a (which in our 
case is mainly the dark energy equation of state w since 
all other parameters are fairly constrained by our priors) , 
one applies eq. |A12j with F' a a, ddi/dpp, Cij and dj all 
computed ignoring magnification bias. To be completely 
explicit, cq. |A12| tells us that the average inferred value 
for parameter p a (e.g. pi = w and so on) is: 

24"= P« + E F ' 9 *l (A14) 



[E ^^[Covfl^CPTC'i) " Cf(zj)) 

+ E M /3">7 Ue -^)" 

7 

where F' 9 is the generalized Fisher matrix (eq. |A10j ) 
(the superscript g is supposed to remind us that magni- 
fication is ignored) i.e. 

and the covariance (with magnification ignored) is: 

[Covf] w = (A16) 

(Cfjzi, Zj ) + 5 l3 /n % )Cf T + Cf(z % )Cf(zj) 
/ sky (2^+l) 

where n% is the mean number of galaxies per unit area in 
redshift bin i. Note that the covariance given above is not 
exact (nor is the implicit assumption that the different i 
modes are uncorrelated), but it should be fairly accurate 
for modes with £ > 10 (recall that we use / s ky = 0.5 
in our worked examples), which dominate the signal-to- 
noise. The same statement can be made about the Fisher 
matrix. 

The fiducial parameter values p Q 's should be chosen to 
be close to the perceived best-fit values — recall that to a 
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modeler who ignores magnification, the perceived best-fit 
differs from the true values. All quantities on the right 
hand side of eq. | A 1 4] . such as F' 9 , Cf T and the co- 
variance, should be evaluated with the parameters set at 
the fiducial values p a 's. The only exception is C™ T (zj) 
which represents the average observed cross-correlation 
and should of course be computed using the true param- 
eter values (e.g. w — — 1 and so on in our thought exper- 
iment in Wip . In principle, the perceived best- fit values 
can differ quite a bit from the true values. We therefore 
apply eq. |A14| iteratively: starting with p a = p^ ue , 
we use eq. |A14| to compute p^ fcrred , and then setting 
Pa = P™ fcrrcd , we reapply eq. |A14| and so on. Generally, 
convergence is achieved with only a few iterations. 
The estimation-bias of interest is p™ ferred — p^ rue . The 



associated (theoretical/model-dependent) error covari- 
ance is given by the inverse of the generalized Fisher 
matrix F' 9 (eq. |A15| ) that ignores magnification, and is 
evaluated at the same p^s as above. 

If one is instead interested in making error forecasts 
when magnification bias is taken into account (as in 
Will) , one should use eq. |A13| with a F' a a that is com- 
puted without ignoring magnification bias i.e. 

where the covariance is given by eq. [22] ■ 
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